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^^ Abstract 

^^ In silicon sensors high densities of electron- hole pairs result in a change of the current pulse shape and spatial distribution 
(S| of the collected charge compared to the situation in presence of low charge carrier densities. This paper presents a detailed 

comparison of numerical simulations with time resolved current measurements on planar silicon sensors using 660 nm 

laser light to create different densities of electron hole pairs. 
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1 1. Introduction 

Silicon sensors are frequently used for the detection of 
radiation. These sensors are built as p-n junctions oper- 
ated under reverse bias. Electron hole pairs are created 
by ionization or direct excitation and induce a current in 
the electrodes when drifting in the electric field until they 
reach a contact. When the created charge carrier densi- 
ties are sufficiently high to modify the electric field in the 
^sensor, significant changes, compared to the situation with 
^low charge carrier densities, are observed (so called plasma 
' effects). These effects have been observed in the detection 
of heavily ionizing particles and with high intensity laser 
light. 

Plasma effects are also expected for experiments at x- 
ray free electron lasers. The studies presented here aim 
at a quantitative understanding of the plasma effects for 
experiments at the European XFEL [T. 

For high charge carrier densities the electrons and holes 
form a so called plasma, which dissolves slowly. The plasma 
boundaries effectively shield its inner region from the ex- 
ternal electric field created by the external bias, thus al- 
tering the induced current pulse and increasing the charge 
collection time [2]. Plasma effects decrease as the elec- 
tric field increases [3]. Using incident ions of different 
masses and energies, the influence of material properties 
on plasma effects has been studied in detail in [4 . 

Electrostatic repulsion effects result in an increased lat- 
eral spread of the collected charge and thus in increased 
charge sharing between pixels, as shown for a-particles in 
[5] . The effects on silicon sensors for the European XFEL 
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were investigated using a focused high intensity laser to 
simulate x-rays in [6l |7] • 

In this work, simulation results are compared to mea- 
sured current pulses showing plasma distortions after illu- 
mination with focused laser light of high intensity. 

2. Experimental setup 

Charge carriers were created with a laser of 660 nm 
wavelength and the time resolved current pulses of the in- 
vestigated diode were read out by a Miteq AM- 1309 wide- 
band amplifier and a Tektronix DPO 7254 2.5 GHz oscil- 
loscope. 

The sample was mounted on a substrate that allowed 
light injection from both sides while providing a stable 
(± 0.1 K, rms) temperature in the range of 240 K to 340 K 
and applying the high voltage to the rear side of the diode. 

The systematic error of the determination of the num- 
ber of generated electron hole pairs was estimated to be 
below 2% by injecting a defined charge into the readout 
system utilizing a defined test capacitance and a voltage 
step function. 

2.1. Equivalent circuit for SPICE simulations 

An equivalent circuit of the setup was used for SPICE 
simulations (see Figure fl]). The equivalent circuit takes de- 
tector capacitance (9.2 pF), signal cable length (3.05 m), 
bandwidth of the amplifier (1 GHz) and imperfections of 
the setup (inductances and capacitances) into account. 
The transfer function of the whole setup has been checked 
and was found to introduce distortions. Thus all simu- 
lations have been convoluted with the transfer function 
(shown in Figure [2|. 
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Figure 1: Equivalent circuit for SPICE simulations. Dominant el- ^g 
ements are the diode capacitance in parallel to the current source 
(marked by the rectangle), the system inductance (30 nH), the lossy 
transmission line (3.05 m RG-174 cable) and the band pass filtering ^o 
properties of the amplifier. si 
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Figure 2: Transfer function derived from the SPICE simulations for 99 
9.2 pF capacitance. The inset shows a zoom of the peak structure of 

the transfer function using the same units on the x- and y-axis. 
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102 

2.2. Laser properties 103 

The laser system [8] emits short and intense hght pulses^°^ 
with a FWHM < 100 ps. The time resolved pulse structure^°^ 
was specified by the manufacturer. For this study laser^°^ 
light with a wavelength of 660 nm (± 2 nm) was used. As^°^ 
the pulse structure depends on the pulse energy, a constant 
energy of 140 pJ was chosen and attenuated with opticahos 
attenuators, which have no effect on the time structure 
of the pulses. 660 nm light has an attenuation length oP^ 
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Figure 3: Geometry of the investigated sample. The diode has a 
thickness of 280 [xm and allows injection of laser light from the junc- 
tion side and opposite to the junction. 



roughly 3 |LLm in silicon at 20° C. 

The laser beam was focused to a spot with a Gaussian 
profile with a = 10 iitm. In air the Rayleigh length (dis- 
tance from focal point to the point where the beam radius 
increases by V^) is approximately 90 |LLm. The high index 
of refraction of silicon (^ 3.6) allows to assume a constant 
lateral beam profile for the entire absorption path. 

2.3. Investigated diode 

The investigated diode was a planar p+nn+ diode fab- 
ricated by CiS [9 . The silicon used is high resistivity n- 
type diffusion oxygenated float zone silicon with <100> 
orientation manufactured by Siltronic [10 . The effective 
doping of the sample was calculated from a capacitance 
measurement as function of voltage and is 8.2x10^^ cm~^. 
The resulting depletion voltage is 49 V with a dark cur- 
rent at the depletion voltage of 0.72 nA. The sample has 
a very low concentration of lattice defects, leading to life 
times in the order of milliseconds, which is many orders of 
magnitude larger than the longest pulse duration recorded 
in this work. Any trapping effects have been considered 
negligible. 

From the measured pad capacitance of 9.23 pF and the 
pad area of (4.95 mm)^ we obtain a thickness of 280 \im. 
using the standard formula for a parallel plate capacitor 
without edge effects. For the measurement of the pad ca- 
pacitance the capacitance of the guard ring to the back- 
plane and the sensor edge is subtracted by the zero ad- 
justment of the capacitance bridge. The remaining effect 
of the guard ring is estimated to be well below 1 %. The 
estimated uncertainty of the diode thickness is ± 2 jam. 

Mechanical measurements of the thickness yield the 
somewhat higher value of 287 \im but include 'dead' layers 
like implantations, passivations and aluminizations. 

The diode has an opening on the p+ side and an alu- 
minum grid on the n+ side to allow light injection. The 
gap between the metalization of diode guard ring is 20 jitm 
wide, the distance between the corresponding implanta- 
tions is 10 |LLm. A sketch of the diode is shown in Figure [3] 

3. Simulations 

For the transport simulations the classical van Roos- 
broeck equations are used [11 . 



Ill 

112 
113 
114 
115 
116 
117 
118 
119 
120 
121 
122 

123 
124 
125 
126 



127 
128 
129 
130 
131 
132 
133 
134 
135 
136 
137 
138 
139 
140 
141 
142 



The lifetime of the plasma cloud is determined by thei44 
emission of carriers into the surrounding volume with lowus 
charge carrier density. A slow movement of the centerue 
of the plasma cloud results from the different probabil-i47 
ity to emit an electron or a hole from different parts ofus 
the plasma, depending on the distance to the closest elec- 
trode. Most of the charge carriers outside of the plasma are 
transported in the depleted volume. The electric field in 
this volume (outside of the plasma) is defined by the bulk 
doping and the external voltage. Dominating influences 
in the simulation are the field induced mobility reduction, 
the initial cloud size and density. 

3.1. Simulated physics 

The fundamental equations describing the problem are 
the Poisson equation and the continuity equations for elec-149 
trons and holes. 150 
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Where e = coCr is the dielectric permittivity, w the elec-^^^ 
trostatic potential, n and p the electron and hole density, ^^g 
C the density of impurities, (j)n,p the quasi- Fermi potential^^g 
for electrons or holes, /in,p the charge carrier mobilities^gQ 
and R the recombination or generation rate. ^^^ 

Potentials are normalized to a constant reference po-^g2 
tential and densities to a reference density. The meaning^g3 
of the variables and parameters used are listed in the ap-^g4 
pendix. ^^^ 

The mobility models used are summarized and dis-^^g 
cussed in [12]. The mobility /i depends on lattice {L),^^^ 
ionized (/) and unionized {N) impurity scattering, carrier-^^g 
carrier scattering {np) and the electric field, where (/), {N)^^^ 
and {np) have the meaning of densities. ^^^ 

Lattice and ionized impurity scattering is modeled us-^^^ 
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Carrier-carrier scattering is implemented using the struc- 
ture of Adler's model [T3l: 
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The total mobility is approximated by using Matthies- 
sen's rule [14 to combine the contributions (they have the 
meaning of scaled inverse macroscopic cross sections E, the 
total cross-section H^o^ is the sum over all S^ of uncorre- 
lated scattering events i). 
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The total mobility is reduced by velocity saturation for 
high values of |V^n,p|- 
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In the simulations the driving force (V0n,p) is used to 
derive the mobility reductiorj^ For the sake of simplicity 
we call this reduction 'field induced'. 

3.2. Numerical methods 

The code follows the theory described in [15]. The 
time integration is based on variable order implicit back- 
ward differentiation formulas (BDF) with time step size 
and order control (see appendix). 

For certain detector applications a very precise num- 
ber of electron hole pairs (relative error << 10~^) has to 
be generated in order to estimate possible charge losses 
in the device. In order to achieve this high precision an 
automated rescaling procedure is used. 

The rescaling procedure uses a comparison of two inte- 
grations. For this, all external generation processes (sources, 
e.g. laser excitation) are parameterized in space and time. 
The first integration is done independent from the simula- 
tion process and results in a high precision source integral 
and the number of electron hole pairs that should be gen- 
erated. 

The second integration is executed during the simu- 
lation and uses the discrete time and space steps of the 
simulation for integration. The integration result is the 
number of generated electron hole pairs in the simulation. 
There is a small difference (<1%) between generated and 
intended number of charge carriers due to the finite step 
sizes and time discretization errors. 

Charge carrier numbers are rescaled to the intended 
numbers, derived values like current are rescaled as well. 
In this way the required high precision generation is as- 
sured. 

The contact currents are evaluated by using test func- 
tions, which approximate the solutions of the related ad- 
joint problem [16]. These techniques are necessary to fulfill 
charge conservation requirements expected in detector de- 
sign applications [17j . 



-•^In the depleted bulk silicon this is identical to the derivation 
using the electric field (Vw). 
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Figure 4: The left graph displays a grid of mesh size h defining the slice of a cylinder. Step sizes are hr = 1, ..., 4:\xm and hz = 0.2, ..., 2\xm. 
Top and bottom of the grid are defined as ohmic contacts, the boundary conditions create a cylindrical symmetry. The right graph displays 
doping profiles as function of depth. The tail of the n+ profile refiects the bulk doping concentration which is marked by the dotted horizontal 



185 4. Simulation parameters 207 

208 

186 A rather simple spatial domain is used to test the ^u.-^^^ 

187 merical methods. A 10° slice of a cylinder of 280 |Ltm height2^Q 

188 and 100 |Ltm radius is discretized by a set of tensor product2^^ 

189 grids that are rotated and split in tetrahedra to obtain di^^^ 

190 slightly anisotropic Delaunay grid. This allows to align the2^3 

191 edges in the main field direction and results in less points 

192 compared to isotropic grids. Thus nested grids, labeled h^^^ 

193 (22990 nodes), h/2 (91339 nodes)and h/4 (364117 nodes) 

194 are possible (see Figure |4|. 

195 Along the main drift path of the charge cloud the grid 

196 is refined (central part of the cylinder). The illuminated 

197 side of the cylinder (top or bottom) is refined as well to 

198 resolve the small initial cloud. The distribution of created 

199 charge carriers is parameterized as 
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with cTx = cTy = 10 |LLm, Zs top or bottom, Xabs = 3 |Ltm the225 
absorption length of the laser light, L{t) the laser pulse226 
shape and Co a constant to ensure that the correct num-227 
ber of electron hole pairs is created. The laser pulses shape 
is shown in the right graph of Figure [5| its peak is very 
narrow (FWHM < 100 ps) and its total length is approx- 
imately 1 ns. 



The p+ doping profile is the result of a one dimensional 
process simulation (compare Figure [4|. It is of short range 
compared to the n+ doping profile. 

The n+ doping profile was determined by spreading 
resistance measurements. The bulk doping concentration 
of 8.2 X lO^^cm"^ is marked by the dotted horizontal line 
in the right graph of Figure |4] 

4.1. Employed mobility parameters 

All results presented here are based on charge carrier 
mobilities reported in literature [12 (summarized in Table 
[2]), but using /3n = 1 (instead of (3^ = 2) for the field depen- 
dent mobility reduction (labeled literature mobility), un- 
less mentioned otherwise. For comparison two other mod- 
els have been used as well. The second model is the same 
as mentioned above but without field dependent mobility 
reduction (labeled constant mobility). The third model 
uses the same models for lattice, ionized and unionized 
impurity scattering as well as the same model for carrier- 
carrier scattering, but the field dependent parameters for 
<100> crystal orientation described in [18] were used (la- 
beled fitted mobility, summarized in Table Isl). 
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Figure 5: The left graph shows the current pulse for different numbers of charge carriers created at the junction side, rescaled to a fixed 
number of charge carriers (1200 electron hole pairs). For cloud sizes with more than 1.2x10^ electron hole pairs pulse distortions due to the 
charge carrier density are visible. The right graph shows pulse shapes for 1200 electron hole pairs starting an the junction (high field side) 
and the influence of the order of the time integration scheme. For comparison the laser pulse structure is also shown. The applied voltage for 
both graphs is 80 V. 



228 5. Simulation results 262 

263 

229 5.1. Test of the numerical methods 

-' 264 

230 To test the numerical methods and the models includ- 

265 

231 ing their parameters a set of simulations with low numbers 

232 of charge carriers was done. The simulated currents were 

233 rescaled to a fixed charge after the simulation process for 

234 comparison and are shown in the left graph of Figure [5] 

235 Typical discretization errors are evaluated by repeat-^^^ 

236 ing computations for one number of charge carriers with^^^ 

237 different spatial and time step sizes and a different time^^° 

238 integration order. Effects of different time integration or-^^^ 

239 ders for junction side creation of charge carriers is shown^^^ 

240 in the right graph of Figure [5J Effects for charge carrier^^^ 

241 creation opposite to the junction are shown in Figure [6)^^^ 

242 The left graph of Figure [6] shows the effects of different ^^^ 

243 time integration schemes, the right graph shows the effect ^^^ 

244 of different grid sizes. 

245 As expected the implicit Euler scheme (i.e. the first or-^^^ 

246 der time integration scheme, see appendix for details) re-^^^ 

247 suits in too much energy dissipation. This can be observed^^° 

248 as additional diffusion and is seen in the pulse shape of the^^^ 

249 currents shown in the right graph of Figure |5] and the left^^^ 

250 graph of Figure [6J 

251 For the implicit Euler scheme the current starts to de-^^^ 

252 crease earlier due to the artificially increased diffusion and^^^ 

253 shows a long tail. 

254 The differences between second and third order time^^^ 

255 integration schemes are small compared to the difference^^^ 

256 between the implicit Euler and the higher order schemes. ^^^ 

257 As a consequence all simulations which were compared^^° 

258 to measurements, do not use the first order scheme in the^^^ 

259 calculations. 

260 The spatial discretization errors are checked by using^^^ 

261 the /i, /i/2, and h/A grids explained above. Compared with^^^ 



the time discretization error of the Euler scheme the space 
discretization error shown in the right graph of Figure [6] is 
small. 

The infiuence of spatial discretization errors is consid- 
ered negligible. 

5.2. Interpretation of the results 

When the charge cloud is created at the junction (i.e. 
high field) side the rise time of the current measures the 
cloud creation and charge carrier separation time. For 
small clouds the current peak corresponds to the situation 
where all charge carriers are drifting. In this case all holes 
are quickly collected at the electrode and the current is 
due to drifting electrons. 

The almost linear reduction of the current is due to the 
lower average velocity at lower electric fields. 

Around the 'knee' the first charge carriers reach the 
contact and the total number of charge carriers in the sen- 
sor drops. The strong decrease of the current represents 
the removal of charge carriers from the bulk volume at 
the contact and thus refiects the shape of the charge cloud 
along the drift direction. This shape is very sensitive to 
the enhanced diffusion of the implicit Euler scheme. 

When the charge cloud is created opposite to the junc- 
tion (i.e. low field) side the current shows a steep increase 
while the cloud it created and separates. In this case all 
electrons are quickly collected at the electrode and the 
current is due to drifting holes. 

The almost linear increase of the current is due to the 
higher drift velocity at higher electric fields. 

Around the peak the charge cloud reaches the contact 
and the number of charge carriers in the bulk material 
drops. Like in the case of illumination from the junction 
side the decrease of the current represents the removal of 
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Figure 6: Comparison of typical pulses of charge carriers created opposite to the junction with different orders in the time integration scheme 
(left) and different grid sizes (right). Both graphs show pulses for 1200 electron hole pairs starting opposite to the junction for an applied 
bias of 80 V. 
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Figure 7: The left graph shows the influence of /3n on pulses for 2.9x10^ electron hole pairs for 80 V and 200 V applied voltage. The right 
graph shows pulse shapes for 1200 electron hole pairs and /3n = 1 as function of the applied voltage. For both graphs the electron hole pairs 
were created at the junction side. 



charge carriers from the sensor volume at the contact andsn 
its shape is very sensitive to the enhanced diffusion of the3i2 
imphcit Euler scheme. sis 

At the apphed voltages neither the drift velocity ofsu 
electrons nor the drift velocity of holes is saturated. Thussis 
a charge cloud of low density traveling in the directionsie 
of a decreasing electrostatic field is compressed, while a 
low density charge cloud traveling in the other direction issi? 
elongated. sis 

The mobility reduction due to large fields has a signif-3^g 
icant infiuence on the pulse shape. A strong infiuence of32Q 
Pn,p is observed as the drift velocity in the entire volume32^ 
is neither in the ohmic regime nor saturated. Figure [7] il-322 
lustrates the infiuence of the model parameter /^^ = 1 or323 
/3n = 2. /3n = 1 results in more pronounced spatial fronts324 
and /3n = 2 causes a lower velocity reduction and thus 



smaller fiight times at high voltages. 

While simulations with f3n = 2 show significantly dif- 
ferent pulses than observed in measurements, simulations 
with /3n = 1 produce pulses which are very similar to the 
measurements, as shown in Figure [9] for junction side illu- 
mination for defocused laser light. 

5.3. Time evolution of the spatial distribution of charge 
carriers 

Figure [S] and the movies in the online version of this 
article show simulations of the time evolution of the hole 
density for 11x10^ electron hole pairs and a bias of 200 V. 
The electron hole pairs are created by laser light of 660 nm 
light focused to 10 |J,m, injected opposite to the p+n junc- 
tion. 



t= 1.8 ns 



t=21.5ns 




■ +5.5 


^ +4 


+2 


■.1 



Figure 8: Evolution of a plasma cloud in space and time. The logarithm to the base of ten of the hole density divided by lO"*^^ cm~^ is 
color-coded for a cut along the y axis. The simulation volume is 100 |J.m high (top to bottom) and 280 |LLm wide (left to right). 11x10^ 
electrons and holes are created on the right side (opposite to the junction) and holes drift to the left. A bias of 200 V is applied. 



325 From the drift-diffusion calculations it is concluded, 

326 that the plasma cloud does not expand, instead charge 

327 carriers are continuously released from the plasma region 

328 and thus form a conductive channel connecting both elec- 

329 trodes. 



330 6. Comparison of simulations and measurements 

331 Electron hole pairs have been created with 660 nm laser 

332 light of different intensities. The measurements for low 

333 intensities allow to verify the simulations, especially the 

334 field dependence of the mobility. For high intensity illumi- 

335 nations plasma effects dominate experimental results and 

336 simulations. 

337 The observed current pulses show small ringing fea- 

338 tures (oscillations), due to the presence of inductances in 

339 the readout circuit. If the gradient on a decreasing slope 

340 is large enough the associated ringing may produce and 

341 undershoot into negative currents. This behavior is mod- 

342 eled by the equivalent circuit shown in Figure [l] The un- 

343 dershoot is neither expected nor found in the simulations 

344 before convolution with the transfer function (Figure [2|. 

345 In the simulations peak densities of np/n'^ ^ 10^^ are 

346 observed, at this density carrier-carrier scattering is still 

347 of small influence. Simulations without mobility reduction^^° 

348 due to carrier-carrier scattering showed the same pulse du-^^^ 

349 rations and are not presented here. Artificially enhanced^^^ 
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Figure 9: Results for junction side illumination using defocused laser 
light to avoid plasma effects. Simulations with literature mobility are 
shown as dashed lines and measurements as solid lines. As simula- 
tions and measurements are very similar the use of /3n = 1 is justified. 



carrier-carrier scattering increases the pulse duration as 
expected, while the measurements show shorter pulses than 
the simulations. 
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Figure 10: Measurements for junction side illumination using focused laser light (solid lines) compared to simulations. Simulations with 
literature mobility are shown as dashed lines, simulations with constant mobility as dotted lines. While acceptable agreement can be observed 
for low densities, high densities show a shorter pulse duration than the simulations. Note the different scales on the time axes. 



353 6.1. Junction side illumination 358 plasma effects. Inspection of the simulations shows peak 

Injection on this side allows to check the transport^^^ values of np/n^^ ^ 10^° in the plasma cloud for the case of 

355 properties of electrons separately from those of holes, as^^° 97xl(j generated charge carriers. 

356 holes reach the close by electrode quickly. ^ei Figures [9] and [lO] show the comparison of simulated and 

Increasing the number of created carriers increases the^^^ measured currents for junction side illumination. 



Simulations and measurements for low density clouds 
are very similar and acceptable agreement between sim- 
ulation and measurement is observed for 1x10^ electron 
hole pairs. Thus it can be concluded that the low density 
transport properties of electrons are reasonably well mod- 
eled. At higher charge carrier densities the current pulse 
shapes deviate from the low density case due to the life- 
time of the plasma cloud. The current is determined by the 
release of charge carriers from the plasma cloud. The de- 
viation between measurements and simulations increases 
with increasing plasma density. 

As seen in Figure [lOJ the current rises slowly at the 
beginning of the pulse for 97x10^ electron hole pairs in 
simulation and measurement. The rise is followed by an 
approximately constant current, until the final decrease 
shows a similar time constant as in the 10x10^ electron 
hole pairs case. 

In the simulations the peak at the beginning is due 
to the removal of the low density periphery of the cloud. 
The remaining high density core has an ellipsoidal shape 
and shrinks with time as well as the maximum density in 
the plasma (see Figure [8|. The barycenter of the plasma 
slowly moves away from the junction. 

Using a constant mobility speeds up the release of charge 
carriers from the plasma and their drift in the rest of the^^^ 
sensor volume. Thus pulses calculated using a const ant^^^ 
mobility are systematically too short, except for the pulse^^° 
obtained for 97x10^ electron hole pairs at 100 V bias,^ 
which is too long by 4% of the pulse length. 

422 

6.2. Illumination opposite to the junction 

This situation allows to study the transport properties 
of the holes. Contrary to electrons, holes move towards the^^^ 
high field region (junction). Figures 



11 



and 12 show the ' 



comparison of simulated and measured currents for this 

427 

case. 

A qualitative agreement between simulation and mea- 
surement is observed for the measurement with defocused^^^ 
laser light and 1x10^ electron hole pairs. However the^^° 
simulated pulses are systematically too long. 

For 11x10^ electron hole pairs the 100 V and 200 V''' 
curves show clear deviations, while at 500 V the deviations 
are at the level observed for 1 x 10^ electron hole pairs. The 
transient behavior between the low density regime and first ^^^ 
plasma effects observed for junction side illumination can^ 
also be observed in this case. 

When 103x10^ electron hole pairs were created, both^ 
the simulated and the measured current pulse become very 
long compared to the current pulses observed for low den-^^° 
sities. While qualitatively the simulations produce pulses^^^ 
with similar durations, the simulated pulses are systemat-^^^ 
ically too long. 

As in the case of illumination on the junction side, 
using a constant mobility speeds up the release of charge^^^ 
carriers from the plasma and their drift in the rest of the^^^ 
sensor volume. Thus pulses calculated using a const ant ^^^ 
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Figure 11: Results for illumination opposite to the junction using 
defocused laser light to avoid plasma effects. Simulations with liter- 
ature mobility are shown as dashed lines and measurements as solid 
lines. 



mobility are systematically too short as well, except for 
the pulse obtained for 11x10^ electron hole pairs at 100 V 
bias (2.4% too long) and 103x10^ electron hole pairs at 
200 V (5.6% too long) and 100 V bias (25.9% too long). 

7. Discussion and Conclusions 

At low charge carrier densities the simulations for the 
collection of electrons (junction side illumination) agree 
well with the measurements. This gives us confidence in 
the simulation as well as in the model used for the electron 
transport. For the collection of holes (opposite side illu- 
mination), the simulated pulses are typically 10% longer 
than the measured ones. This may indicate a problem 
with the simulation of the charge cloud separation at low 
fields or with the model used for the hole transport. In 
[18 the experimental data has been used to determine the 
parameters of electron and hole transport. 

High charge carrier densities result in significant dis- 
tortions of the pulse shapes, both in data as well as in 
simulations. For the highest charge densities investigated 
(?^ 10^ electron hole pairs) the measured pulse lengths in- 
crease by up to a factor 5 for electrons and a factor 10 
for holes. The results are charge collection times up to 
250 ns for a detector operated well above depletion! The 
simulations qualitatively describe the data and provide an 
understanding of the pulse shape for high charge densi- 
ties: An initial rise due to the collection of charges from 
the periphery of the charge cloud is followed by an approx- 
imately constant current due to the release of charges from 
a shrinking plasma cloud. 

A number of attempts have been made to understand 
the discrepancies. 
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Figure 12: Measurements for injection opposite to the junction using focused laser light (solid lines) compared to simulations. Simulations 
with literature mobility are shown as dashed lines, simulations with constant mobility as dotted lines. Measurements systematically show a 
shorter pulse duration than the simulations. Note the different scales on the time axes. 



To estimate the effect of having an initial charge car-454 
rier distribution which is different from the hght profile, 455 
simulations with different widths of the initial charge car-456 
rier distribution and absorption lengths have been done457 
for 1x10^ and 11x10^ electron hole pairs created opposite458 



to the junction. The results for 200 V applied bias are 



presented in the upper graphs of Figure [13] and [Mj show- 
ing that an increase in width by 50 % or an increase in 
absorption length by 100 % is not sufficient to produce a 
pulse which is as short as the measurement. 
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NeM [10^ 



1 
1 
1 

10 
10 
10 
97 
97 
97 

1 

1 

1 

11 
11 
11 
103 
103 
103 



Uuas [V] 



500 
200 
80 
500 
200 
100 



500 
200 
100 
500 
200 
100 
500 
200 
100 
500 
200 
100 
500 
200 
100 
500 
200 
100 



duration [nsl 



0/7.0/6.8 
0/8.5/8.4 
0/13.0/13.5 
0/9.0/10.0 
0/13.5/15.0 
0/23.0/27.0 



5.0/7.0/7.0 

6.5/9.0/9.0 

10.0/12.0/12.3 

5.0/7.0/7.0 

8.0/11.0/12.0 

15.0/17.0/20.0 

8/14/18 

24/28/38 

52/50/75 

7.0/9.0/10.0 

13/15/17 

31/31/38 

9/12/14 

27/28/35 

87/85/100 

25/30/35 

95/90/120 

340/270/390 



difference 



/ -2.9% 
/ -1.2% 
/ 3.8% 
/ 11.1% 
/ 11.1% 
/ 17.4% 



-28.6% / 0.0% 
-27.8% / 0.0% 
-16.7% / 2.5% 
-28.6% / 0.0% 
-27.3% / 9.1% 
-11.8% / 17.6% 
-42.9% / 28.6% 
-14.3% / 35.7% 

4.0% / 50.0% 
-22.2% / 11.1% 
-13.3% / 13.3% 

0.0% / 22.6% 

-25% / 16.6% 
-3.6% / 25.0% 

2.4% / 17.6% 
-16.7% / 16.7% 

5.6% / 33.3% 
25.9% / 44.4% 



Table 1: Pulse durations for all measurements (tmeas) presented in this work compared to their simulations with the literature mobility 
model (tut) and the constant mobility model (tconst)- The column labeled Ne,h lists the number of created electron hole pairs; a '-' indicates 
the low density measurement and simulation. The column labeled Uhias lists the applied bias voltage. The column labeled duration lists 
tconst /tmeas /tut- The column labeled difference lists {tconst - tmeas) /tmeas and {tm - tmeas) /tmeas- A '0' indicates that no corresponding 
simulation was performed. The measurements have been performed with 660 nm light focused to a Gaussian spot with criaser = 3 |Ltm on a 
p+nn+ diode with a thickness of 280 fxm and an effective doping of 8.2x10-'^-'^ cm~^. 
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Using literature data of the mobility, the low density 483 



limit for holes (Figure 11) could not be reproduced, indi-484 
eating the need for changes in the mobility parameteriza-485 
tions. Using the fitted mobility parameterization for holes486 
the current pulse for 1x10^ electron hole pairs could be487 



well reproduced (Figure 13, bottom). 



Although certainly not realistic, we have made Simula- 489 
tions assuming that the mobility is independent of the elec-49o 
trie field (constant mobility model). The pulses obtained49i 
are systematically shorter than the measured ones (except492 
for the highest intensity at low bias voltage). However493 
the mobility model from the literature and the constant494 
mobility model give the interval in which the measured495 
pulse durations are found. Table [l] summarizes the re-496 
suits. The results assuming the constant mobility modeW 
are also presented in Figs. [lO] and [T2J and in the upper498 
graphs of Figs. [13] and [l4j They show, that the release of499 
charge carriers from the plasma cloud can be significantlysoo 
accelerated, when the diffusion parameter, which is relatedsoi 
in our simulation to the mobility by the Einstein relation502 
{Dn^p = l-in^p^BT/q)^ does not decrease with increasing503 
field. Combined with the increased drift velocity of the504 
charge carriers this results in almost all cases in simulated 
current pulses, which are shorter than the measured ones. 



This observation suggests that the Einstein relation be- 
tween mobility and diffusion, which is valid only in the 
equilibrium case, has to be modified at high fields and/or 
at high charge carrier densities. 

The possible effects of Fermi-Dirac statistics (see [22] . 
[23] ) were estimated and the expected infiuence on the 
pulse length should be studied, especially because the Ein- 
stein relation {D = /j^ksT/q) would be replaced by a den- 
sity dependent one. Hence the interplay of drift and diffu- 
sion will be modified. 

The van Roosbroeck equations used in the simulations 
describe the situation in the diode for low density charge 
clouds well, however in situations with high densities gra- 
dients may violate assumptions made in the derivation of 
the drift-diffusion approximation (compare [12] for a short 
discussion). 

Optical photons of 660 nm wavelength (1.87 eV) pro- 
duce so called hot charge carriers, as the band gap of silicon 
is approximately 1.12 eV. The simulation assumes charge 
carriers in thermal equilibrium with the crystal lattice, 
which is justified, as the thermalization of the hot carriers 
is usually very fast compared to the pulse duration. 
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Figure 13: Effects of different parameters on the simulated current pulses for 1x10^ electron hole pairs at 200 V applied voltage. The upper 
graph shows the influence of different mobility models. The lower graph shows the influence of different distributions of the initial charge 
cloud. Although lower density clouds dissolve faster, the effects of different attenuation lengths are small compared to the effects of a constant 
mobility model. 

505 8. Summary sos with emphasis on the impact on sensor performance of de- 

509 tectors for experiments at the European XFEL. 
A simulation program was developed to model the trans-^^ The numerical stability and apphcability for sensor de- 

507 port of charge carriers for high densities in silicon sensors 
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Figure 14: Effects of different parameters on the simulated current pulses for 11x10^ electron hole pairs at 200 V applied voltage. The upper 
graph shows the influence of different mobility models. The lower graph shows the influence of different distributions of the initial charge 
cloud. Although lower density clouds dissolve faster, the effects of different initial distributions are small compared to the effects of a constant 
mobility model. 

511 sign purposes has been demonstrated. su fects cannot be described by using the mobihty and diffu- 

512 As a result of the comparison of measurements andsis sion models in literature. It is shown that, except for the 

513 simulations it is concluded that the observed plasma ef-sie highest intensity, two different sets of mobility models can 
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Parameter 


Value 




er 


11.67 




/^° 


1448 cmVVs 




M° 


495 cmVVs 




OLn 


2.33 




Oi-p 


2.23 




I 


rii 




c-/ 


3x10^6 cm-3 




c-/ 


4x10^^ cm-3 




^n 


350 




Op 


81 




N 


Ei\Ci\ 




^Adll 


1.04x1021/ T \3/2 
cmVs ^}300K^ 


543 


f^Adl2 


7.45x10^^/ T \2 
cm2 V300K/ 


544 


Vsatr, 


1.1 X 10^ cm/s 


545 


Vsatp 


0.95 X 10^ cm/s 


546 


Pn 


1 or 2 


547 


Pp 


1 


548 



Table 2: Parameters used in the mobility parameterization. Param- 
eters are quotes from il2j, /3n has been set to unity. ^^° 



517 be used to simulate pulses which are either systematically 553 

518 longer or systematically shorter than the measurementS554 

519 and thus allows to estimate the minimum and the maxi-555 

520 mum of the pulse duration. 

521 In spite of the discussed discrepancies the simulatioussg 

522 program is a valuable tool for the design and optimization 

523 of sensors and readout electronics for the European XFEL. 

524 Combining the effects of Fermi-Dirac statistics in the 

525 diffusion process, variations in the parameters of the initial 

526 distribution of charge carriers and an optimized mobility 

527 model it may be possible to provide simulations which re- 

528 produce the measurements reasonably well for all intensi- 

529 ties with a single, unified set of parameters. 
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536 10. Appendix A 

537 10.1. Simulation parameters and derived values 

538 All symbols and their meaning are listed in Table [4j All 

539 potentials in this table are normalized to a constant refer-557 

540 ence potential {Ut = ksTo/q ^ 1/40V at room tempera-sss 

541 ture), all densities to a reference density (n^g/ = lO^^cm"^ ^ 



Hi). 



560 
561 
562 
563 
564 



Parameter 


Value 


VsatpiT) 


A'7A cm-' rp -2.619 
474 y^ Irel 

0.94 X 10^^ Trer°-^^^ 

1.181 Trel°-^'^ 


VsatjT) 

lin(T) 


1440 ""{;' Trel'^"^^ 
1.054 X 10^2^ Trer°-^°^ 

0.992 T'e/°-"^ 



Table 3: Parameters from 18 used for the parameterization labeled 
'fitted mobility'. The abbreviation T^ei = Tqok ^^^ been used to 
improve legibility. 



10.2. Assumptions on the simulation domain, grids, and 
space discretization 

The simulation space is defined on a bounded, polyhe- 
dral domain Q C M^ with a boundary dft = Td UTn- ^d 
denotes the Dirichlet part of the boundary. It is closed, 
has a positive measure and describes ohmic contacts. The 
Neumann part T^ describes insulating boundary parts or 
symmetries. The solutions {w{t)^n{t)^p{t)) are defined in 
S X Q with the time interval 5* = {O^tend) and {w{0) = 
wq^ n(0) = no,p(0) = po) the initial values. The domain is 
discretized by a boundary conforming, tetrahedral Delau- 
nay mesh and the usual Scharfetter-Gummel discretiza- 
tion is used (compare [15], [E]). 

10.3. Time integration 

For the implicit Euler scheme dissipativity can be pro- 
ven. It guarantees positive solutions on any Delaunay grid 
and for any time step. This cannot be expected for higher 
order schemes in general. Hence the related family of back- 
ward differentiation methods (BDF) is used to reduce dis- 
sipation while the order control handles the observed os- 
cillations in space and time. Let 



y'{t) = f{t,y{t)), y{to) 



Vo 



(13) 



denote a nonlinear first order initial value problem with 
the initial value hq. The BDF formulas for order m = 
1, 2, 3 (m = 1 implicit Euler) at time step y{tk) = Vk with 
constant step size r are: 



m = 1 : yu+i - Vk = ^/(^fe+i, ^fe+i) 



(14) 



3 1 

m = 2: -Vk+i - 2yk + -yk-i = rf{tk+i,yk+i) (15) 



2^ 

-^Vk+i - ^yk -^ 7:yk-i 
2 



-yk-2 = rf{tk+i,yk+i) 
(16) 



For 771 = 1,2 the formulas are A-stable (stable for the 
linear test problem and eigenvalues in the left half plane, 
the cases m = 3,4 are still stable for a sector of 88°, 72° 
around the negative real axes), see [24 and literature cited 
therein. The BDF formulas can be modified for variable 
step size. Time step size control is based on predictor- 
corrector differences of functionals, like free energy, dis- 
sipation rate, one selected contact current, sources, and 
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Parameter 




Meaning 


t 






time 


X 






position in space (vector) 


eo 






dielectric constant 


e = eoCr 






dielectric permittivity 


kB 






Boltzmann constant 


Q 






elementary charge 


mo 






electron mass 


^n,p 






effective mass of electron or hole 


asor 






Bohr radius 


h 






reduced Planck constant 


T 






lattice temperature 


To = 293.15i^ 




reference temperature 


w 






electrostatic potential 


(t)n=W- 


log n/ui 




quasi-Fermi potential for electrons 


(l)p =w^\ogp/ni 




quasi-Fermi potential for holes 


n = n^e^" 


-Ckr^ 




electron density 


p = me'^p' 


-w 




hole density 


C 






density of impurities 


rn = rp = 


lO-^s 




Shockley-Read-Hall recombination lifetime 


Po = no = 


rii 




intrinsic charge carrier density 


^ - r.po4 


1 f^2 _ 

-Tpno-\-Tnp-\-Tpn^ i 


-np) 


recombination / generation rate 


Mn,p > 






carrier mobilities 


J-^n,p ^ f^n 


,pkBT/q 




Einstein relation 



Table 4: Parameters and derived values used in the simulations. 
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580 
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their deviations from a polynomial predictor and the num-593 
ber of Newton steps needed to solve the nonlinear equa-594 
tions. The variable time integration order is controlled by 595 
testing different order predictors against the computed so-596 
lution at the present time step. The order of integration597 
of the next time step is defined by the predictor with min-598 
imal error. Time step rejections are based on the local599 
truncation errors and the number of Newton steps. A re-eoo 
jection is combined with order reduction. The maximumeoi 
order can be specified. The identical technique is used to602 
apply continuation methods (the continuation parametereos 
{t) replaces time (t)) with respect to the boundary val-604 
ues or model parameters, because it is often impossibleeos 
to reach large applied voltages directly from the uniquelyeoe 
defined equilibrium, t = corresponds to a state with607 
known solution, t = 1 to the wanted state with an un-eos 
known solution. The only difference, with respect to time609 
integration, is taking truncation errors not into account, eio 
because a sequence of stationary solutions is constructed. en 
The predictor is chosen such that the positivity of the den-612 
sity variables is preserved. 613 

Each time step requires the precise solution of one non- 
linear system of equations by Newton's method. Charge 
conservation is lost accordingly due to the errors intro- 
duced in the solution of the nonlinear equations. eis 



of all dependencies, including the derivatives of parame- 
ters, together. Hence, common expressions are reused in 
the function and the Jacobian. The solution of the nonlin- 
ear equations is controlled by the L'^-norm of all potential 
updates, hence it is precise for low densities, too. 

The linear systems are solved by a combination of di- 
rect [20 , iterative methods [21 and primary and secondary 
preconditioners. The primary preconditioner avoids the 
factorization of the complete Jacobian by approximating 
the main dependencies by scalar equations on the whole 
domain. The secondary preconditioners solve local sys- 
tems to include the missing couplings (avalanche, small 
time steps, strong recombination etc.) for all variables at 
a given grid point or very small subsets of the grid. 

This allows to avoid the factorization of the complete 
Jacobian, hence to reduce memory requirements and the 
operations needed while one is still getting quadratic con- 
vergence of the Newton process. 

The main algorithms defining the limits of the appli- 
cations are grid generation and the solution of the linear 
systems. 
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